Optimized experimental settings for the best detection of quantum nonlocality 
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Nonlocality lies at the core of quantum mechanics from both a fundamental and applicative point 
of view. It is typically revealed by a Bell test, that is by violation of a Bell inequality, whose success 
depends both on the state of the system and on parameters linked to experimental settings. This 
leads to find, given the state, optimized parameters for a successful test. Here we provide, for a quite 
general class of quantum states, the explicit expressions of these optimized parameters and point 
out that, for a continuous change of the state, the corresponding suitable experimental settings may 
unexpectedly vary discontinuously. We finally show in a paradigmatic open quantum system that 
this abrupt "jump" of the experimental settings may even occur during the time evolution of the 
system. These jumps must be taken into account in order not to compromise the correct detection 
of nonlocality in the system. 
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Violation of a Bell inequality indicates, in two-qubit 
(spin-l/2-like) systems, the presence of quantum cor- 
relations nonrcproducible by any classical local model 
(nonlocal correlations) [H H]- A violation is identified 
by values of the Bell function B, consisting of a combi- 
nation of correlations averages, larger than a maximum 
classical threshold i?ciass- For a given system the Bell 
function does depend both on the state of the system 
and on parameters linked to spin-like measurement di- 
rections that correspondingly fix the experimental set- 
tings. Many experimental tests have proven the occur- 
rence of Bell inequality violations in different contexts as 
entangled-polarized photons [21 SI [S] , photon-atom sys- 
tems [6], atom-atom systems [7] and recently also su- 
perconducting Josephson phase qubits . Bell inequali- 
ties have also been extended to characterize the presence 
of nonlocality in multiqubit systems (Mermin inequali- 
ties) j9j [To] . Nonlocal correlations also play a crucial role 
in quantum information science, as in the realization of 
device-independent and security-proof quantum key dis- 
tribution protocols pn [T^ . 

For pure states of any degree of entanglement it is pos- 
sible to show that, by appropriately choosing the Bell 
parameters, violation of Bell inequality always occurs 
|13j . This is not the case for mixed states which may 
present entanglement without any corresponding direct 
violation of a Bell inequality [TJ]. Indeed, mixed states 
take place as a consequence of unavoidable environmen- 
tal noise which causes decoherence and a nonlocal cor- 
relations decay during the evolution [TSJ [TB]. In order 
to test at best the presence of nonlocality, it is essential 
to adopt optimized Bell parameters (OBPs) so that the 
maximum value of the Bell function -Bmax is measured. 
These OBPs in turn determine the appropriate experi- 



mental settings which allow to perform the optimal Bell 
test in order to state that, if iJmax < Bciass, quantum 
nonlocality is surely nondetectable directly. Because of 
the dependence of the OBPs on the state of the system, 
these are expected to change in time. If one wants to 
exploit nonlocal features of the system at time t, it is 
required to detect at best the presence of quantum non- 
local correlations at that time. Therefore, one needs to 
know how these OBPs vary in time so to correspondingly 
adjust the experimental settings and measure i?max(i)- 
A procedure to obtain the OBPs and the correspond- 
ing maximum value of the Bell function for an arbitrary 
two-qubit state is known |17j . However, the explicit ex- 
pressions of these OBPs have never been reported before 
but for some very limited class of pure quantum states 

[I1[IH]. 

Here we firstly give the explicit expressions of the 
optimized Bell parameters for a quite general class 
of pure and mixed quantum states showing moreover 
that, changing the quantum state in a continuous way, 
the OBPs may unexpectedly undergo a discontinuous 
change. These expressions are given for a two-spin- 1/2 
system where the parameters are the angles fixing direc- 
tions in real space. Nevertheless, these parameters can be 
translated into the suitable experimental settings for any 
two-qubit system realization in different physical frame- 
works. We finally consider the inescapable environmen- 
tal effects on the evolution of the two-qubit system and 
their influence on the Bell test. In a paradigmatic open 
quantum system we show that the abrupt discontinuous 
change of the OBPs may occur during the time evolution 
and this must be considered in order not to compromise 
the correct detection of nonlocality in the system. 



I. QUANTUM NONLOCALITY 



'Elect ronic address: [Io franco@fisica. u"nipa.it| 
tURL: http: //www.f isica.unipa. it/~compagno 



In the following we use the Clauser-Horne-Shimony- 
Holt (CHSH) form of Bell inequality [1 [S], which is 



2 



the most suitable for an experimental test of nonlocal- 
ity in bipartite quantum systems composed by spin-1 /2 
objects. Let the operator Os — OsiOs,4's) be a spin 
observable with eigenvalues ±1 associated to the spin 
5" = 1,2 defined by Os — O5 • erg, where O5 = 
(sin cos (/)s, sin 6*5 sin 0s, cos 6*5) is the unit vector in- 
dicating an arbitrary direction in the spin space and 
as = {a^ ,ay ,a^) is the Pauli matrices vector. The 
CHSH-Bell inequality associated to a two-spin-1 /2 state 
p is then B{p) < 2 [21 [17l US], where B{p) is the Bell 
function defined as 

B{p) = \(0,02) + (0,0',) + (0[02) - {0[0',)l (1) 

where (0102) = Tr{/50iC'2} is the correlation function 
of observables d, O2 and O'^ = Os{0's, <Ps)- given the 
state p, it is possible to find a set of angles {6i,9'i,92, ^2} 
and (/)2, (^2} such that the CHSH-Bell inequality 

is violated {B{p) > 2) the correlations are nonlocal and 
cannot be reproduced by any classical local model. 

A necessary and sufficient condition for violating 
the CHSH-Bell inequality by an arbitrary two-spin-1/2 
mixed state is known [T7j. Once introduced the 3x3 
real matrix Tp, defined by the matrix elements tmn = 
Tr(/5(T„ (E) (Trn), whcrc the cr's are the Pauli matrices. 



one must build the matrix U„ 



a symmet- 



ric matrix which can be therefore diagonalized {Tj in- 
dicates the transposed matrix of Tp). The maximum 
of the Bell function B{p) results then to be given by 
-Bmax(/5) = 2^maXj>fc{Mj -I- Wfc}, where j, fc = 1, 2, 3 and 
Uj's are the three eigenvalues of the matrix Up. A proce- 
dure to obtain the angles giving the maximum value of 
the Bell function is also provided ^ . 



II. MAXIMUM OF BELL FUNCTION 

In this section we give the maximum value of the 
Bell function for the class of states whose density ma- 
trix px, in the standard computational basis B = {|1) = 
111), |2) = |10), |3) EE |01), |4) = |00)}, has a X structure 
of the kind 



PX 
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(2) 



This class of states is sufficiently general to include the 
two-qubit states most considered both theoretically and 
experimentally, like Bell states (pure two-qubit maxi- 
mally entangled states) and Werner states (mixture of 
Bell states with white noise) [10» .ISj |20j • Moreover, a X 
structure density matrix arises in a wide variety of phys- 
ical situations [3 [HI [22 [23l [25 . A further remarkable 
aspect of the X states is that, under various kinds of dy- 
namics, the initial X structure is maintained during the 
evolution 12011211. 



Using the criterion previously described to obtain 
Smax, the three eigenvalues m's, in terms of the density 
matrix elements, are found to be 

Ui = 4(|/9i4| -I- |p23|)^, U2 = {pil + P44 " P22 - Pss)^, 
U3=4(|pl4|-|p23|)^ (3) 

as also already reported [261 . Being ui always larger than 
M3 , the maximum of Bell function for X states results to 
be 



III. 



-Bmax(/5) = 2Jui + maXj=2,3{"i}- 



OPTIMIZED BELL PARAMETERS 



(4) 



Here we give the OBPs, that is the angles 6's and 0's 
maximizing the Bell function B{p) of Eq. for the X 
states class of Eq. Q . Following a known procedure [T7] 
we find the explicit OBPs, which result in two different 
sets of angles in correspondence of the two regions U2 > 
M3 (region 1) and U3 > U2 (region 2). 

Set 1: U2 > U3. The OBPs in this case are given by 



7r/2, 



0, 



- - sign(pii 



P44 - P22 - P33) arctan 



t/2 = 7r-y2, 

(jji = - [arg(pi4) + arg(p23)] /2, 0'i = 0, 

(1)2 = 02 = [arg(p23) - arg(pi4)] /2, (5) 

where the the function sign (a;) is +1 if x > and —1 if 
X < 0. 

Set 2: U3 > M2- The OBPs in this case are 



Wi = 02 = ^2 = 7r/2, 
- [arg(pi4) -l-arg(p23)] /2, 
01 -Hsign(|p23| - |pi4|)7r/2, 
arg(p23)-arg(pi4) _^^^^^^ /^^ 

2 V 
arg(p23)-arg(M _ /^^ 
2 V "1 



(6) 



Let us observe that the angles of the two sets are not 
equal when ^2 = '"3. This means that, when this condi- 
tion is satisfied, two different OBPs and so experimental 
settings do give the same value of -Bmax- Moreover, con- 
sidering the OBPs as function of the quantum state, to 
a slight continuous change of the quantum state crossing 
the point U2 — U3, it corresponds an unexpected finite 
"jump" of the suitable experimental settings. 

It is worth to emphasize the importance to have an an- 
alytical result for the two-qubit configuration here con- 
sidered. In fact, if an analytical criterion to find OBPs 
was absent, numerical procedures would be necessary. 
This is the case, for example, in systems of more than 
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two qubits where Mermin inequaHties are considered to 
test nonlocahty [10]. On the basis on the two-qubit re- 
sult, one may expect that a finite jump of the OBPs take 
place even for the multiqubit configuration. However, 
if this jump was present, starting from one of the OBPs 
sets, numerical procedures could give rise to serious prob- 
lems in identifying the point of this possible OBPs jump 
because of the finite distance between the two sets of 
parameters. 

It is clear that the angles of Eqs. Q and (|6|, which de- 
termine the orientations of the measurement apparatus 
in real space giving i?niax for two spin-1/2 objects, do cor- 
respond to different experimental settings for a general 
two-qubit system. In cavity QED, for example, qubits 
are represented by either two- level Rydberg atoms or cav- 
ities having one or zero photons; Ramsey zones followed 
by field ionization detectors are then the devices suitable 
for measuring the qubit along an arbitrary pseudospin di- 
rection, whose angles 9, (j> are controlled by the classical 
microwave field amplitude and the atom-field interaction 
time inside the Ramsey zone [27l [28l [29l [30] . In the con- 
text of traveling polarized photons, where the CHSH-Bell 
inequality was originally applied |31j . the spin levels are 
usually coded in the horizontal (H) and vertical {V) po- 
larizations of the photon, the angles cf)'s are usually set 
at zero and the angles 9's are related to the angles 9p's 
of photon polarizers by 9p = 9/2 [2J ITS]; the measure- 
ment is finally performed by photodetectors. In solid- 
state physics with Josephson phase qubits, coded by the 
two first levels of a well potential, pseudospin directions 
are selected by rotating the qubit states using microwave 
pulses before performing measurements by a pulse of cur- 
rent along the z axis of the Bloch sphere (basis states of 
the qubit) [5]. 



IV. OPTIMIZED BELL PARAMETERS IN A 
DYNAMICAL CONTEXT 

So far we have found the optimized parameters (OBPs) 
that give the best condition for the Bell test on a given 
two-qubit X state. From an experimental point of view, 
one firstly needs to know the density matrix of the state 
on which to measure Bmax and confirm if there is a vio- 
lation of the CHSH-Bell inequality. This can be done by 
quantum state tomography Therefore, in presence 
of a time-evolving quantum state, one may establish the 
form of the state with its density matrix elements at the 
time of the measurement. The experimental settings of 
the devices for the Bell test can be then appropriately 
adjusted according to our indications. In this section we 
give a dynamical example where the OBPs not only are 
crucial to reveal nonlocality, but they also present an un- 
expected behavior. 

Consider a system of two identical independent qubits, 
namely 1 and 2, each embedded in a zero-temperature 
bosonic environment. Each qubit has ground and excited 
states |0), |1) separated by a transition frequency luq. For 



our study we need the two-qubit reduced density matrix 
at the time t. This can be accomplished by exploiting 
a procedure based on the knowledge of the single-qubit 
dynamics |25j . The dynamics of the single part "qubit- 
environment" is described, in the case of bosonic reservoir 
under rotating wave and dipole approximations, by the 
Hamiltonian f35] 

H = fi.wo(7+(7_ + ^fi ^ujkblbk + [gk^+bk + gla^bl^ , 

k 

(7) 

where cr+ = |1)(0|, (j_ = |0)(1| are the qubit raising and 
lowering operators, 6^, bk are the photon creation and 
annihilation operators and gk is the coupling constant 
of the mode k with frequency LOk- It is found that, for 
a zero-temperature environment, the single-qubit evolu- 
tion under the Hamiltonian of Eq. Q can be described 
in terms of a single function q(t), whose square modulus 
|g(t)P can be linked to the single-qubit excited state pop- 
ulation. In the basis {|1), |0)}, the single-qubit reduced 
density matrix p^{t) can be then written as f 32j 

V piimit) Po^o(o)+pfi(o)(i-k(t)n; 

(8) 

The fact that the density matrix p^{t) can be expressed in 
terms of a single function q{t) is general, the information 
on the environment spectral density and coupling con- 
stants being contained in the explicit time dependence 
of the function q(t) itself. The two-qubit density matrix 
p(t) is then obtained [25j, its elements depending only on 
their initial values and on the function q{t), regardless 
of the environment structure. Hereafter, for the sake of 
simphcity, we shall indicate |q(t)P = 

We choose as two-qubit initial state the extended 
Werner- like state defined as I^Uj 

p*(0) = r|<i>)(<i>| + (l-r)/4/4, (9) 

where r is the purity parameter, /4 the 4x4 identity 
matrix and |<I>) — a\01) + /3e*''|10) is the one-excitation 
Bell- like state with a,/3 positive real and + — 1. 
/5*(0) reduces to a Werner state for a = /S = 1/^/2 and 
(5 = 0, TT, that is when its pure part becomes a Bell state. 
For r = the state above becomes a totally mixed state, 
while for r = 1 it reduces to the Bell-like (pure) state |<I>). 
This state has an X form which is maintained during the 
dynamics due to the spin-boson Hamiltonian of Eq. ([7| . 
We can analyze the dependence of the dynamics of -Bmax 
on the initial mixedness and degree of entanglement. 

Starting from the two-qubit state defined above, the 
expressions of Ui, U2 and of Eq. ([3]) are found to be 
in terms of jqp as 

u, = U3^ 4aWkl', ^^2 = [1 - 2|g|2 + (1 - r)\q\^]\ 

The condition U2 = determines a closed surface in the 
3D space (a^, r, Igp) given in Fig.n] with two values jgip. 
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FIG. 1: Surface of the quantum state condition U2 = «3 in 
the (a-^,r, space. Given a^,r two values of \q\^ (klP and 
|(jr2p) satisfy U2 = U3. The points outside the surface represent the 
region 1 where 112 > «3, while those inside the surface represent 
the region 2 where U3 > U2. 



|(72p for a fixed couple a^, r. The points enclosed by this 
surface represent the region 2 (u^ > U2), while the points 
outside represent the region 1 (u2 > u^). Dynamically 
|(7p changes with time and Fig.fllclearly shows that, given 
and r, when |gp crosses the surface at time t, the 
OBPs immediately before and after this time are given 
respectively by the two different sets 1 or 2. This in 
turn gives place to a discontinuous change with time of 
the experimental settings giving the maximum of Bell 
function. 

In order to examine this surprising aspect, we plot in 
Fig. [2] both Bi = 2y/ui + U2 and B2 = 2y/ui + U3 as a 
function of \q\'^ for — 0.3 and r = 1. In this case 



the initial state of Eq. ([9]) is a pure non-maximally en- 
tangled state and -Bmax = max{i3i, B2}. In the plot one 
moves from right to left with time going from 1 to 
0) and i?inax is given by Bi (branch 1; solid curve) before 
|(7(A)p and by B2 (branch 2; dashed curve) after it. The 
region of violation of CHSH-Bell inequality (-Bmax > 2) 
is zoomed in the inset. At value |q(A)p there is the pas- 
sage from branch 1 to branch 2 and the corresponding 
abrupt discontinuous "jump" of the OBPs (see part b 
of Fig. 2]) which pass from set 1 to set 2. Moreover, for 



values 



9(C)|^ 



< 



< 



|(7(B)|^, using the parameters of 



set 1 (giving the maximal violation when \q\ > |g(A)| ) 
one finds no violation of Bell inequality while using those 
of set 2 one finds a violation. 

The environment determines the explicit time evolu- 
tion of |gp = |g(i)p. For Markovian environments |(?(t)P 
decays exponentially [32] and the values |g(A)p and 
|(7(B)p are typically reached after a short time. Differ- 
ently for structured environments inhibiting the popula- 
tion decay, such as photonic crystals, |q'(t)P decay can be 
slow [33] and the crossing value |(7(A)p could be reached 
after a rather long time. Another case is that of non- 
Markovian environments where collapses and revivals of 
|(7(t)P may be present, such as in high-Q cavities [25j . 
and the surface of Fig. [l] could be crossed several times 
during the evolution. 

Therefore, when one needs to detect the presence of 




|q-(A)|^ 



|q+(A)|^ 



FIG. 2: Dynamical behavior of Bi and B2 with abrupt 
"jump" of the OBPs. a, Plot of Bi 



and B2 



2^u\ + U2 (solid curve) 
2\/"l + "3 (dashed curve) as a function of \q\-^, for 



a = 0.3 and r = 1. The inset evidences the violation region 
where B\,B2 > 2. The value |g(A)p indicates the change from 
branch 1 to branch 2 (going from right to left), while between the 
values |g(B)|2 and |(?(C)P only B2 = ^max is larger than two. b, 
The spin-1/2 measurement vectors determined by the OBPs are 
illustrated (for (5 = 0) in correspondence of |g(A)p. The abrupt 
jump of the OBPs is well evidenced while passing from |g+(A)p to 
|(;_(A)p. Immediately before and after the crossing point |(ir(A)p, 
the spin measurement vectors lie on the x-z plane and x-y plane, 
respectively. 



nonlocal correlations without ambiguity, the striking be- 
havior of the OBPs here shown highlights the importance 
of our analysis to fix the most appropriate experimental 
settings to reveal nonlocality. 



V. CONCLUSIONS 

In this paper we have found the optimized parameters 
determining the experimental settings suitable to reveal 
at best the presence of nonlocality in a bipartite spin-1 /2- 
like system. We have reported the explicit dependence on 
the state of these optimized Bell parameters (OBPs) for 
a general class of quantum states. We have found that 
this class of states can be divided in two parts each of 
them corresponding to a different set of OBPs. We have 
also found that even if the state of the system changes 
continuously, as when it evolves with time, a surprising 
abrupt jump of these OBPs may occur in proximity of 
some states. In particular, we have shown that these 
jumps of the experimental settings that maximize the 
Bell function occur during the evolution of a paradig- 
matic two-qubit open system. This phenomenon must 
therefore be taken into account in any experimental test 
aiming at obtaining the best confirmation of the presence 



of nonlocality in a quantum system. 



[1] J. S. Bell, Physics 1, 195 (1964). 

[2] J. F. Clauscr, M. A. Home, A. Shimony, and R. A. Holt, 
Phys. Rev. Lett. 23, 880 (1969). 

[3] P. G. Kwiat, S. Barraza-Lopez, A. Stefanov, and 
N. Gisin, Nature 409, 1014 (2001). 

[4] S. Groblacher et al., Nature 446, 871 (2007). 

[5] A. Fedrizzi et al.. Nature Phys. 5, 389 (2009). 

[6] D. L. Moehring, M. J. Madsen, B. B. Bhnov, and C. Mon- 
roe, Phys. Rev. Lett. 93, 090410 (2004). 

[7] M. A. Rowe et al.. Nature 409, 791 (2001). 

[8] M. Ansmann et al.. Nature 461, 504 (2009). 

[9] N. D. Mermin, Phys. Rev. Lett. 65, 1838 (1990). 
[10] R. Horodecki, P. Horodecki, M. Horodecki, and 

K. Horodecki, Rev. Mod. Phys. 81, 865 (2009). 
[11] A. Acin, N. Gisin, and L. Masanes, Phys. Rev. Lett. 97, 
120405 (2006). 

[12] N. Gisin and R. Thew, Nature Photon. 1, 165 (2007). 

[13] N. Gisin, Phys. Lett. A 154, 201 (1991). 

[14] R. F. Werner, Phys. Rev. A 40, 4277 (1989). 

[15] B. BcUomo, R. Lo Franco, and G. Compagno, Phys. Rev. 

A 78, 062309 (2008). 
[16] B. Bellomo, R. Lo Franco, and G. Compagno, Adv. Sci. 

Lett. 2, 459 (2009). 
[17] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 

Lett. A 200, 340 (1995). 
[18] M. A. Nielsen and I. L. Chuang, Quantum Computation 

and Quantum Information (Cambridge University Press, 

2000). 

[19] J. F. Clauser and M. A. Home, Phys. Rev. D 10, 526 
(1974). 

[20] B. Bellomo, R. Lo Franco, and G. Compagno, Phys. Rev. 

A 77, 032342 (2008). 
[21] E. Hagley et al., Phys. Rev. Lett. 79, 1 (1997). 
[22] S. Bosc, I. Fuentes-Guridi, P. L. Knight, and V. Vedral, 

Phys. Rev. Lett. 87, 050401 (2001). 
[23] J. S. Pratt, Phys. Rev. Lett. 93, 237205 (2004). 
[24] J. Wang, H. Batelaan, J. Podany, and A. F. Starace, J. 

Phys. B 39, 4343 (2006). 
[25] B. Bellomo, R. Lo Franco, and G. Compagno, Phys. Rev. 

Lett. 99, 160502 (2007). 
[26] L. Derkacz and L. Jakobczyk, Phys. Rev. A 72, 042321 

(2005). 

[27] C. C. Gerry, Phys. Rev. A 53, 4583 (1996). 

[28] R. Lo Franco, G. Compagno, A. Messina, and A. Napoli, 

Phys. Rev. A 72, 053806 (2005). 
[29] R. Lo Franco, G. Compagno, A. Messina, and A. Napoli, 

Phys. Rev. A 76, 011804(R) (2007). 
[30] S. Haroche and J. M. Raimond, Exploring the Quantum: 

Atoms, Cavities, and Photons (Oxford University Press, 

USA, Oxford, New York, 2006). 
[31] A. Aspect, P. Grangier, and G. Roger, Phys. Rev. Lett. 

49, 91 (1982). 

[32] H.-P. Breucr and F. Petruccione, The Theory of Open 
Quantum Systems (Oxford University Press, Oxford, 
New York, 2002). 

[33] B. Bellomo, R. Lo Franco, S. Maniscalco, and G. Com- 
pagno, Phys. Rev. A 78, 060302(R) (2008). 



